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Abstract 

We study Christoffel and Darboux transforms of discrete isother- 
mic nets in 4-dimensional Euclidean space: definitions and basic prop- 
erties are derived. Analogies with the smooth case are discussed and a 
definition for discrete Ribaucour congruences is given. Surfaces of con- 
stant mean curvature are special among all isothermic surfaces: they 
can be characterized by the fact that their parallel constant mean cur- 
vature surfaces are Christoffel and Darboux transforms at the same 
time. This characterization is used to define discrete nets of constant 
mean curvature. Basic properties of discrete nets of constant mean 
curvature are derived. 

1 Introduction 



Stimulated by the integrable system approach to isothermic surfaces (cf. [|J 
or H) — or, better: to Darboux pairs of isothermic surfaces^ - - "discrete 
isothermic surfaces" or "nets" have been introduced by Alexander Bobenko 
and Ulrich Pinkall ||. For the development of this theory it seems hav- 
ing been crucial not to use the standard calculus for Mobius geometry but 

* Partially supported by the Alexander von Humboldt Stiftung and by the Forschungsin- 
stitut fur Mathematik, ETH Zurich 

1 Transformations seem to play an important role in the relation of surface theory and 
integrable system theory. 
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a quaternionic approach — which is well developed in case of Euclidean 
ambient space. In || the Christoffel transform (or "dual") |7j for discrete 
isothermic nets in Euclidean space — the Christoffel transform of an isother- 
mic surface is closely related to the Euclidean geometry of the ambient space 
- is developed. The Christoffel transform is an important tool in the the- 
ory of isothermic surfaces in Euclidean space: it can be used to characterize 
minimal surfaces (cf. 0, ||) and surfaces of constant mean curvature (cf. 
[|I~2]j ) between all isothermic surfaces — for discrete nets of constant mean 
curvature we will later use this characterization as a definition. 

Recently, a quaternionic calculus for Mobius differential geometry [|I(J was 



developed. This approach turned out to be very well adapted to the theory 
of (Darboux pairs of) isothermic surfaces [TT| - based on these results a 



quaternionic description for the Darboux transform of isothermic surfaces in 
Euclidean 4-space J? 4 = H was given WQ: any Darboux transform of an 
isothermic surface in f? 4 can be obtained as the solution of a Riccati type 
partial differential equation. This equation can easily be discretized which 
provides us with a discrete version of the Darboux transform. Note that it 
seems to be more natural to work in 4-dimensional ambient space rather than 
in the codimension 1 setting^: for example the structure of our Riccati type 
equation becomes more clear in 4-dimensional ambient space. 

After a comprehensive discussion of the cross ratio in Euclidean 4-space 
we recall the definition of a discrete isothermic net and prove the basic facts 
on the Christoffel transform in 4-dimensional ambient space. Then, we give a 
definition of the Darboux transform of discrete isothermic nets by discretiz- 
ing the Riccati type partial differential equation which defines the Darboux 
transforms of an isothermic surface in the smooth case. In the smooth set- 
ting a Darboux transform of an isothermic surface is obtained as the second 
envelope of a suitable 2-parameter family of 2-spheres. In the discrete set- 
ting there naturally appears a "2-parameter family" of 2-spheres, too — this 
suggests a definition of the "discrete envelopes" of a "discrete Ribaucour 
congruence" which we briefly discuss. Following the definitions we prove 
Bianchi's permutability theorems [TJ for multiple Darboux transforms and 
for the Darboux and the Christoffel transform in the discrete case. In the 



2 Apart from the fact that calculations can be done much easier in 3- or 4-dimensional 
ambient space it seems to be clear that most considerations on the Darboux transform 
of discrete isothermic nets can be done in arbitrary dimensions — ironically, there might 
occur problems with the Christoffel transform. 
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final section we discuss discrete nets of constant mean curvature: (smooth) 
cmc surfaces can be characterized by the fact that their parallel cmc surface 
is a Christoffel and a Darboux transform of the original surface at the same 
time. We use this as a definition and discuss Darboux transforms of constant 
mean curvature of discrete nets of constant mean curvature. 

To illustrate the effect of the Darboux transformation on discrete isother- 
mic nets we have included a couple of pictures. Comparing these pictures 
with those in |12|] (which were calculated using the "smooth theory") the 
reader will observe considerable similarities — indicating another time the 
close relation with the "smooth theory". A discrete isothermic net in Eu- 
clidean 3-space allows oo 4 Darboux transforms into discrete isothermic nets 
in 3-space. Figure [I] shows a Darboux transform of a discrete isothermic 
net on the Clifford torus which is closed in one direction. The reader will 
recognize a typical behaviour of Darboux transforms: along one curvature 
line "bubbles" are "added" to the original net while ,in the other direction, 
the transform approaches the original net asymptotically. From this typical 
behaviour we might expect that Darboux transforms of (smooth or discrete) 
isothermic nets on a torus will never become doubly periodic. 

Figure ^ shows a periodic constant mean curvature Darboux transform of 
a discrete isothermic net on the cylinder — the correspoding smooth Darboux 



transform is a Backhand transform of the cylinder, at the same time [0 . The 
(also periodic) isothermic net (of constant mean curvature) shown in figure 
|lOl is obtained by applying a second Darboux transformation. In the last 
section of the present paper we will show that discrete isothermic nets of 
constant mean curvature H ^ allow oo 3 Darboux transforms of constant 
mean curvature H. A similar theorem holds for discrete minimal netsQ (cf. 



— a minimal Darboux transform of an isothermic net on the catenoid 



is shown in figure [?[ Other (obviously non minimal) Darboux transforms 
of the same catenoid net with "bubbles added" along different directions are 
displayed in figures [| and — these surfaces might help to "see" the "typical 
behaviour" of the minimal Darboux transform in figure [7| . . . 

3 Discrete minimal nets are defined in ||. As the proof for the theorem on cmc Darboux 
transforms of cmc nets relies on the permutability theorem for multiple Darboux trans- 
forms, the proof for that on minimal Darboux transforms of minimal nets depends on the 
permutability theorem for Christoffel and Darboux transforms — when we characterize 
minimal isothermic nets by the fact that their Christoffel transforms are isothermic nets 
on a 2-sphere. 



3 




Figure 1: A typical Darboux transform of the Clifford Torus 



Since the cross ratio plays a key role in our investigations on discrete 
isothermic surfaces and the Darboux transform it seems useful to discuss 
properties of 

2 The Cross Ratio 

of four points in Euclidean 4-space. In H Alexander Bobenko and Ulrich 
Pinkall introduce the (complex) cross ratio 

DV(Q 1; Q 2 , Q 3 , Qa) = ReQ(Qi,Q 2 , Q 3 , Q 4 ) ± « • |ImQ(Q 1; Q 2 , Q 3 , Q 4 )| 

where Q(Qi,Q 2 ,Q 3 ,Qa) = (Qi ~ Q 2 )(Q 2 - Qs^iQs - Q 4 )(Q 4 - Qi)' 1 for 
four points Qi ) Q2 ) Qz)Qi G R 3 — ImU in Euclidean 3-space. If S C R 3 is 
a 2-sphere containing these four points this is exactly the cross ratio of the 
four points interpreted as complex numbers on S as the Riemann sphere^ - 

4 Note that this cross ratio is slightly different from the classical cross ratio: the classical 
one may be obtained by interchanging Q2 and Q3 in the above formula. This reflects the 
fact that in discrete surface theory it is more useful to consider the cross ratio as an 
invariant of a surface patch rather than an invariant of two point pairs. 
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the ambiguity of the sign of the imaginary part in the cross ratio corresponds 
to the possible orientations of the sphere. Note that the imaginary part of 
the cross ratio vanishes exactly when the four points are concircularQ In this 
case the four points do not determine a unique 2-sphere but a whole pencil 
of 2-spheres. 

This construction holds perfectly for four points Qi, % = 1 . . .4, in 4- 
dimensional Euclidean space R 4 = H: 

Definition (Cross Ratio): The complex number 

BV{Q 1 ,Q 2 ,Q 3 ,Q i ) := ReQ(Qi, Q 2 , Q 3 , Qa) + i- |ImQ(Q 1; Q 2 , Q 3 , Q 4 )|, ( ] 
Q(Qi,Q2,Q3,Q±) := (Qi-Q2)(Q2-Q 3 )-\Q3-Qa)(Qa-Qi)- 1 1 ' 

is called the CROSS RATIO of four points Qi,Q 2 ,Q 3 ,Q4 <E R 4 = H in Eu- 
clidean 4-spacef\ — or of the quadrilateral (Qi, Q 2 , Q 3 , Q 4 ) C R 4 , respectively. 

It is a straightforward calculation to express the cross ratio in terms of 
the distances of the four points[} 

[Z?2^34 + ^14^23 ~~ ^13^24] + H ~ ^((hj)i,j=l,-A) 

W(Q 1 ,Q 2 ,Q 3 ,Q i ) = -g-z^ (2) 

zi 14 i 23 

where Uj := \Qi — Qj\- Note that the vertices of a planar quadrilateral are 
concircular if and only if the product of the diagonal lengthes equals the sum 
of the products of the lengthes of opposite sides|. Consequently, the cross 
ratio of concircular points can be calculated from the lengthes of two pairs 
of opposite edges. 

5 Or, they arc colinear. — But since we are doing Mobius geometry we do not distinguish 
circles and straight lines. 

6 Note that for four points in the complex plane C C H (which already carries an 
orientation) we may simply use Q itself as the cross ratio. 

7 At this point we would like to thank Matthias Ziircher for helpful discussions. 

8 This is reflected by the fact that 

- det((J?.)i,j=x,...,4) = (a + b + c)(-a + b + c)(a-b + c)(a + b-c) 

where a = I12I34, b = I13I24 and c = luhs- Note that all factors are non negative which 
can be interpreted as triangle inequalities for the triangle with lengthes a, b and c. With 
this interpretation, the determinant in the cross ratio is 16 times the squared area of this 
triangle (Heron's formula). 
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If we consider J? 4 = H = {(Q, 1) | Q G -H"} as (a subset of) the conformal 
4-sphere HP 1 = {V ■ H \ V G H 2 } the squared norms in this equation 
are replaced by a biquadratic function \Qi — Qj\ 2 = T>((Qi, 1), (Qj,l)), i.e. 
V{V\,W) = V(V,WX) = \\\ 2 V(V,W) for V, W G H 2 and Ae£ Thus 
the cross ratio of four points is independent of the choice of homogeneous 



coordinates of the points and consequently it is a conformal invariant |T0 . 

As another consequence of the above equation (|2|) we derive the following 
identities for the cross ratio — note that the determinant in (0) is invariant 
under permutations of the four points: 

DV(Q 4 ,Qt,Q 2 ,Q 3 ) = m{ Qi ^ QzQA) , 
DV(Q 1) g 3 ,Q 2 ,Q4) = 1-W(Q 1 ,Q 2 ,Q 3 ,Q 4 ). 

The complex conjugation in these equations arises since we defined the cross 
ratio to have always a non negative imaginary part. With the help of these 
two identities we easily obtain a second set of identities, 

DV(Q 1 ,Q 2 ,Q 3 ,Q 4 ) = DV{Q 3 ,Q A ,Q 1 ,Q 2 ) 

= DV(Q 2 ,Q 1 ,Q 4 ,g 3 ) (4) 
= DV(Q 4 ,Q 3 ,Q2,Qi). 

The complete set of the twenty four identities for the cross ratio can be 
derived using these two sets of identites — they are symbolized in figure 0. 

As a preparation we first give a planar version of our hexahedron lemma 
— which will take a central role throughout this paper: 

Lemma: Given a quadrilateral (xi,x 2 ,x 3 ,X4) in the complex plane and a 
complex number A G C there is a unique quadrilateral (z%, z 2 , z 3 , Z4) to each 
initial point Z\ G C such that the cross ratios^ satisfy 

Q{Z!, z 2 , z 3 , Zi) = Q{x 1 ,x 2 ,x 3 ,x 4 ) =: 11 

QOi, 22,^2,^1) = QO3, Zi, x 4 ,x 3 ) = fi\ (5) 

Q,(Z 2 ,Z 3 ,X 3 ,X 2 ) = Q(Z4, Z\,X\,X±) = A. 



This lemma may be proven by a straightforward algebraic computation: 
first we solve the equations Q(zi, z 2 , x 2 , Xi) = fi\, Q(z 2 , z 3 , x 3 , x 2 ) = A and 

9 Here, we use Q as the cross ratio of four points in the complex plane (cf. footnote @). 
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Figure 2: Identities for the complex cross ratio 
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Q(zs, ^4) ^4) #3) — A*A succesively for z 2 , -23 and Z4. This defines us a unique 
quadrilateral (zi, Z2, za,z&) in the complex plane. What is left to do is to 
verify the remaining two cross ratios — since the calculations become fairly 
long we prefered to let the computer algebra system "Mathematica" do the 
work: 

QDVSolve [a_ , Ql_ , Q2_ , Q3_] := Block [{X}, 
X = (Q2-Q3) (Q1-Q2) A (-1) a; 
(1+X) A (-1) (X Ql + Q3) 

]; 

(* cross ratio of the initial quadrilateral: *) 

x4 = QDVSolve[ m, xl, x2, x3 ]; 

(* we compute the second quadrilateral: *) 

z2 = QDVSolve[ m 1 , x2 , xl , zl ] ; 

z3 = QDVSolve [ 1 1, x3, x2, z2 ]; 

z4 = QDVSolve[ ml, x4, x3, z3 ] ; 

(* we verify the remaining two cross ratios: *) 

Factor [ zl - QDVSolve [ 1, xl, x4, z4 ] ] 

Factor [ z4 - QDVSolve [ m, zl, z2, z3 ] ] 

This completes the proof of the above lemma — note, that we could have 
omitted the last calculation since the cross ratio of the constructed quadri- 
lateral can also be derived directly from figure |||. 

This lemma can not be directly generalized to quadrilaterals in space since 
a given (complex) cross ratio does no longer determine the fourth vertex of a 
quadrilateral in space. But, in our applications we will only be interested in 
real cross ratios — and in this case the cross ratio does determine the fourth 
vertex. Thus, if we start with an initial quadrilateral (Xi, X 2 , X3, X4) with 
concircular vertices and one vertex Z\ of the other vertex we can construct 
succesively the other vertices Z 2 , Z 3 and Z 4 on the circles given by the 
point triples {Xi,X 2 , Z{\, {X 2 , X 3 , Z 2 } and {X 3 , X4, Z 3 }. Since these circles 
are always given by three points lying on the 2-sphere S containing the 
circle through the vertices of the first quadrilateral and the initial vertex Z\ 
the whole construction takes place on this 2-sphere. Consequently, we can 
prescribe one (real) cross ratio as above: 

The Hexahedron lemma: Given a quadrilateral (Xi, X 2 , X%, X4) in Eu- 
clidean 4-space with concircular vertices and a real number A £ R there is a 



8 



z 4 




Figure 3: The hexahedron lemma (planar version) 

unique quadrilateral (Zi, Z 2 , Z 3 , Z4) to each initial point Z\ e H such that 

DV(Z 1 ,Z 2 ,Z 3 ,Z 4 ) = DV(X 1 ,X 2 ,X 3 ,X 4 ) =: n 
DV(Z 1 ,Z 2 ,X 2 ,X 1 ) = BV(Z 3 ,Z 4 ,X,,X 3 ) = //A (6) 
DV(Z 2 ,Z 3 ,X 3 ,X 2 ) = W{Z^Z U X U X±) = A. 

Moreover, the vertices of the hexahedron (Xi, X 2 , X 3 , X 4 ; Z l5 Z 2 , Z 3 , Z±) lie 
on a 2-sphere S C H. 

Now, we are prepared to start our investigations on 

3 Discrete isothermic nets 

Away from umbilics isothermic surfaces in Euclidean 3-space can be defined 
by the existence of conformal curvature line parameters; the results from 
1 1] and [|l^] suggest that this is a good definition in case of 4-dimensional 
ambient space, too. In the codimension 2 case this definition requires the 
normal bundle of the surface to be flat — otherwise, it makes no sense to 
speak of "curvature line coordinates" . 
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Figure 4: A Darboux transform of the Catenoid 



Following 0] a surface is parametrized by curvature lines if and only if 
the parameter lines divide the surface in infinitesimal (planar) rectangles and 
the surface is isothermic if and only if it is divided into infinitesimal squares 
by its lines of curvature (cf. ||). We may reformulate these criteria in a 
Mobius invariant flavour using the cross ratio: 

Lemma: A (smooth) surface is parametrized by curvature lines if and only 
if the parameter lines devide the surface in infinitesimal patches with negative 
real cross ratios and 

the surface is isothermic if and only if its curvature lines devide the surface 
in infinitesimal harmonic patches, i.e. with cross ratios -1. 

These criteria may motivate the following definition of discrete isothermic 
nets — note that we prefer the notion of a "net" since this term refers to a 
parametrized surface rather than to a surface in space as a purely geometric 
object: a "discrete isothermic net" is the analog of an isothermic surface 
given through conformal curvature line parameters. 
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Definition (Discrete Curvature line Net, Discrete Isothermic Net): 

A map F : T — > H, T c I 2 , is called a discrete curvature line net if 
the vertices of all elementary quadrilaterals are concircular and the quadri- 
laterals are embedded: 

) < 0; (7) 

it is called an isothermic net if all elementary quadrilaterals are harmonic: 

(8) 



It is possible to define more general "isothermic nets" or "isothermic sur- 
faces" as a discrete version of isothermic surfaces given in arbitrary curvature 
line coordinates (see H). Even though some formulas have to be modified 
in this approach (see for example the definition of the "dual surface" in ||) 
all the facts on the Christoffel and Darboux transforms of isothermic nets we 
are going to develop seem to hold in the more general setting. However, for 
this presentation we prefer the simpler and more enlightening approach we 
just introduced. 

As in the codimension 1 case, discrete isothermic nets have a Christoffel 
transform (dual). And, in some sense, they can be characterized by the 
existence of a Christoffel transform: 



Theorem and Definition (Christoffel transform): If a discrete net 
F : T — ► H is isothermic then 

A G R \ {0} ; defines another discrete isothermic net F c : T — > H , a 
Christoffel transform (or dual,) of F. Note that a Christoffel trans- 
form F c of an isothermic net F is determined up to a scaling and its position 
in space. 

If, on the other hand, (fjj defines a second discrete net — i.e. ^) is integrable 
— then both nets are isothermic, or parallelogram nets. 
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This theorem has an analog for smooth surfaces: given a parametrized 
surface F : U C M 2 — > H the integrability condition for the //-valued 1- 
form F~ x dx — F~ l dy is equivalent to F xy = and F conformal — i.e. F 
is a conformal curvature line parametrization of an isothermic surface — or 
F xy = and hence F(x, y) = Fi(x) + F 2 (y) is a translation surface — which 
is the smooth analog of a parallelogram net. 

During the proof of this theorem let us denote by a := F m+l n — F m>n , 

b ■ = fm+l,n+l •^m+l,n) C ■ = ^m,n+l ^m+l,n+l and .= F m ^ n F m ^ n +i the 

edges of an elementary quadrilateral of the net F. Clearly, we have 

= a + b + c + d (10) 

— the closing condition for the net F. The closing condition ( "integrability" ) 
for the dual net F c reads 



= i-i+l-I 

a_ b _ c a 

a b I c 

— \a\ 2 \b\ 2 ^ |c| 2 \d\ 2 



a ^ b C | _5 d_ (11) 

UI2 IM2 ~T UI2 IJI2 ■ 



As a first consequence this shows that a necessary condition for the net F c to 
close up is that the elementary quadrilaterals of F be planar. If we assume 
the quadrilateral of F to be planar we can rewrite ( |TT] ) as 

= a (I_i + l_i) c 

= -[1 + ab~ 1 cd- 1 }d- [1 + ad- 1 cb- 1 ]b (12) 
= [1 + ab^cd" 1 ] ■ [a + c] 



where we used bed = deb for the last equality: see fll^p . 

If F is an isothermic net its elementary quadrilaterals are planar with 
cross ratio akr x cd~ x = —1. Consequently, the net F c closes up. The only 
thing left is to show that elementary quadrilaterals of F c have cross ratio -1. 
But, 

-1 = ab^cd" 1 = d- x cb- x a = \\b\d}- 1 (13) 

which shows that F c is an isothermic net, too. This proves the first part of 
our theorem. 



To attack the second part we note that from the above calculation (12) it 
follows that two cases can occur if is integrable: the quadrilateral under 
investigation has cross ratio — 1 = ab~ l cd~ l or it is a parallelogram. Assum- 
ing that the quadrilaterals of F do not change their type we conclude that 
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F and (from the first part) F c are isothermic nets or both nets are parallel- 
ogram nets. This completes the proof of the above theorem — provided we 
happen to prove the following 



Lemma: Two vectors a, b G R 4 = H are linearly dependent if and only if 

ah = ba; (14) 
three vectors a,b,c G H are linearly dependent if and only if 

abc = cba; (15) 
four vectors a,b,c,d G M are linearly dependent if and only if 

abed — deba = deba — abed. (16) 

The first statement is immediately clear since this simply means that 
a t G R. To prove the second statement we calculate 

abc — cba = 2[det(ai, b±, ci), — (ao • b\ x c\ + b ■ c\ x a\ + Co • a\ x b\)\ 

where Xq = Hex and X\ = Imx. A careful analysis of this equation provides 
the second statement. And third we have 

Re[abcd — deba] = 2 det(a, b, c, d) 

which completes the proofQ 

To motivate our ansatz for the Darboux transform of discrete isothermic 
nets let us shortly recall some facts on 

4 The Darboux transform 

of smooth isothermic surfaces. A sphere congruence (a 2-parameter family 
of spheres) S : M 2 — > {spheres and planes in R 3 } is said to be "enveloped" 
by a surface F : M 2 — ► R? if at each point the surface has first order contact 

10 At this point we would like to thank Ekkehard Tjaden for helpful discussions. 
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to the corresponding sphere[]]: F{j>) G S(p) and d p F(T p M) = Tp( p )S. If a 
sphere congruence S has two envelopes F and F — which is, generically, the 
case — then it establishes a point to point correspondance between its two 
envelopes. From the works of Darboux and Blaschke |2| we know that if 
this correspondance preserves curvature lines — the sphere congruence is a 
"Ribaucour sphere congruence" — and it is conformal, too, then genericallyQ 
both envelopes are isothermic surfaces. In this case the two envelopes are said 
to form a "Darboux pair" , one surface is said to be a "Darboux transform" 
of the other. These definitions can be generalized to codimension 2 surfaces 
in 4-dimensional space by introducing congruences of 2-spheres — which are 



not hyper spheres any more [11 



In [|12| we derived a Riccati type partial differential equation 

dF = (F- F)dF c (F - F) (17) 

for Darboux transforms F of an isothermic surface F where F c is a Christoffel 
transform of F. Note that in order to obtain every Darboux transform F of 



F as a solution of (15), it is crucial not to fix the scaling of the Christoffel 
transform: rescalings of F c make F run through the associated family of 
Darboux transforms. Considering this Riccati equation as an initial value 
problem we see that an isothermic surface allows oo 5 Darboux transforms 
- or, if we are interested in surfaces in 3-dimensional space, it allows oo 4 
Darboux transforms. 

This Riccati type equation (|r7|) can now be easily discretized to obtain a 
system of difference equations 

^(-fm+l,n Fm,n^) (^'rn.n Fm,n){^F m +l,n Fm,n) ( ^'m + l.,n ^m+l,n) 

•^(-^m,w+l ^m,n) {Fm,n Fm,n) (-^m,n+l -^m,n) (^m,n+l -^m,n+l) 

where we replaced the Christoffel transform F c of F according to (|9|). How- 
ever, this ansatz is not unique: we also could have interchanged the roles of 
the edges connecting F and F by replacing the previous equations by 

^(-fm+l,n Fm,n) (^ m - 1 ,n ^tn+l,n)(^m+l,ii Fm,n) ( f'r nv F mn } 

•^(-^m,n+l ^*m,7i) (^m,n+l F'm,n-\-l) (-^m,n+l Fm,n) (-^m,n -^m,n) 



11 For simplicity reasons we do not give exact definitions — we just want to give the ideas 
of the used terms. For a comprehensive discussion the reader is refered to the classical 
book of W. Blaschke [§. 

12 If the sphere congruence is, in a certain sense, "full" , i.e. if it is not a congruence of 
planes in some space of constant curvature. 
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Figure 5: Another Darboux transform of the Catenoid 



- or we even could have used a mean value of the two. But, rewriting 
these equations with the cross ratio we see from (f|) that both ansatzes are 
equivalent. Thus, without loss of generality we may use our first ansatz to 
obtain the following 



Theorem and Definition (Darboux Transform): If F : T — > H is an 

isothermic net then the Riccati type system 



^ DV(-^m,ri) F mnj -fm+l,n> -^m+l,n) 



(18) 



defines an isothermic net F : V — > H. Any solution F of ( fiq j is called a 
Darboux transform of F. 

This theorem is an easy consequence of the hexahedron lemma (page |9|). 
Note that from the identities (|4]) for the cross ratio we see that the equations 
(|18D are symmetric in F and F — consequently, we also may refer to the 
pair (F, F) of isothermic nets Darboux pair. 
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As a second consequence from the hexahedron lemma we conclude that 
for a Darboux pair (F, F) the vertices of the hexahedron 

(-^m,ni ^m+l,ti) -F«.+l,n+l) -F?i,n+1] Fn,nj ^m+l,ti) ^m+l,n+l i ^m,n+l) 

lie on a sphere S( mjn y — these spheres do naturally live on the "dual lattice" 

T* := {((m,ra),(m+l,n),(m+l,n + l),(m,ra+l))| , . 

(ffl,n),(m+l,n),(ra+l,n+l),(m,rj + l)er} 1 J 

of T which consists of the elementary quadrilaterals. Obviously, an (interior) 
point pair (F m , n ,F m , n ) can be obtained as the intersection of four spheres 
S( m , n )*, ^(m-i.n)*, ^(m-i.n-i)* and ^(m.n-i)* ■ These facts suggest that the 
Darboux pair (F, F) "envelopes" a discrete Ribaucour congruence: 



Definition (Envelopes of a Discrete Ribaucour Congruence): Two 

discrete curvature line nets F, F : T — > I? 4 are sazd to envelope a dis- 
crete RIBAUCOUR CONGRUENCE S : T* — > {2-spheres and 2-planes in R A } 
if (interior) point pairs (F m ^ n , F m ^ n ) lie on four consecutive spheres S( m ,n)*> 

S(m-l,n)*, ^(m-l.n-l)* an d S^n-l)* ■ 

Note that we require the two envelopes to be curvature line nets — and 
consequently the sphere congruence to be Ribaucour. If we would omit this 
assumption the patches of one envelope could determine the spheres of the en- 
veloped congruence uniquely since the four vertices of an elementary quadri- 
lateral would not be concircular in general. But this would contradict the 
fact that in the smooth case there is always a pointwise 1-parameter freedom 
for an enveloped sphere congruence when one envelope is given. 

The previous discussions indicate a strong similarity to the smooth case — 
to complete this picture we prove two fundamental "permutability theorems" 
which hold in the smooth case: 



Theorem: If Fi,F 2 :T — > H are two Darboux transforms of an isothermic 
net F : F — > H with parameters Ai and X 2 , respectively, then there exists an 
isothermic net F = F i2 = F 2i :Y^H which is a X 2 -Darboux transform of 
Fi and a Xi-Darboux transform of F 2 at the same time. The nets F, F\, F 2 
and F have constant cross ratio 

£ = DV(F,F 2 ,F,A). (20) 
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F 



Figure 6: Bianchi's permutability theorem 



This theorem is a discrete analog of Bianchi's permutability theorem M. 
Again, it can be obtained as a consequence of the hexahedron lemma: let 
us pick one edge of an elementery quadrilateral of F and the corresponding 
edges of F\ and F 2 , as indicated in figure || We denote the vertices of these 
edges by X,Y, Xi,Yy and X 2 ,Y 2 . Since Fx and F2 are Darboux transforms 
of F we have 

TN(X,X X ,Y U Y) = ±\ 1} 
DV(X,X 2 ,Y 2 ,Y) = ±A 2 . 

Then, we can find two points X, Y such that 

dv(x,1 2 ,x,1 1 ) = DV(y,y 2 ,y,y 1 ) = £ 

and, by the hexahedron lemma, 

DV(A > 1 ,x,y,y 1 ) = ±a 2 , 

DV(X 2 ,X,Y,Y 2 ) = ±Al 

Repeating this construction for all four edges — note that by construction the 
resulting edges close up to form a quadrilateral — we obtain an elementary 
quadrilateral of an isothermic net JP which is a Ai-Darboux transform of F 2 
and a A2-Darboux transform of F\, as desired. This completes the proof of 
our first permutability theorem. 
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Since the four surfaces in the permutability theorem have constant cross 
ratio (|20| ) it becomes clear that there is an even fancier version of the per- 
mutability theorem which involves not just four but eight surfaces (cf. |T[): 
if we start with three Darboux transforms Fi, F 2 and F 3 of an isothermic net 
F we can construct three isothermic nets F\ 2 , F 23 and F 3 % via 

£ = DV(F,F 2 ,F 12 ,F 1 ), 

% = DV(F,F 3 ,F 23 ,F 2 ), (21) 
£ = DV(F,F 1 ,F 31 ,F 3 ). 

Now, we may apply the construction a second time to obtain an eighth net 
F which satisfies 

£ = BV(F 3 ,F 3U F,F 23 ), 

£ = FN(F 1 ,F 12 ,F,F 31 ), (22) 
£ = DV(F 2l F 23 ,F,F 12 ). 

It is a consequence of the hexahedron lemma that the net F exists. This 
seems to be worth formulating a 

Corollary: If Fi : V —>■ H , i = 1,2,3, are Darboux transforms of an 
isothermic net F : F — > H with parameters \, then there exist four isother- 
mic nets Fij, F : T — ► H , ij = 12, 23, 31, such that each F^ is a \-Darboux 
transform of Fj and a \j -Darboux transform of Fi and such that F is a X%- 
Darboux transform of F 23 , a X 2 -Darboux transform of F 3 \ and a X 3 -Darboux 
transform of F\ 2 , at the same time. 

The second "permutability theorem" is a discrete version of the compat- 
ibility theorem for the Darboux and Christoffel transform of an isothermic 
surface pf (cf. 0): 

Theorem: If F, F : T — > H form a Darboux pair, then their Christoffel 
transforms F c , F c = F c : F — >• H form — if correctly scaled and positioned 
— a Darboux pair, too. 

To prove this theorem we imitate the proof we gave in the smooth case 
[|T2]]: first, let us consider the difference functions G := F — F : F —> II and 
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G c := F c — F c : T — > H. Since F is a Darboux transform of F we have 

•^(-^m+ljTj Fm,n) (-^m,n -^m,n) {-^m+l,n -^m,n) (-^m+l,n ^m+l,n) 

Subtracting the first equation from the second one yields 

1 1 \ / 1 1 



A 



Fm+l,n -^m+l,n F m ^ n F m ^ n J \ F m +\^ n F m n -^m+l,n Fm,r 

and consequently 



G m +l,n G m ^ n 

if F c and F c are scaled as in (H). In a similar way we obtain the same 
difference equation for the other direction 

/~tc _ /--yc ^ ^ 

Thus, we can position F c and F c in R 4 such that Gm,n = ^J„. As in (|T3l), 
we can now calculate the cross ratios 

TW( F c F c F c F c ) — ^- — A 

^ v V- 1 m,n' 1 m,n' 1 m+l,n' x m+l,n) \ 

DVfF c F c F c F c ) = — A 

v V m,ni m,m m+l,m m+l,n) 

showing that _F C and F c form indeed a Darboux pair. This completes the 
proof. 

Troughout the rest of the paper we will consider 



5 Discrete nets of constant mean curvature 

As is well known, (smooth) surfaces of constant mean curvature (cmc) in 
3-dimensional Euclidean space are isothermic surfaces, i.e. - away from 
umbilics — they allow conformal curvature line parametrizations. In fact, 
we know that surfaces of constant mean curvature (H ^ 0) are very special 
examples of isothermic surfaces in the context of Darboux and Christoffel 
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transforms: the (correctly scaled and positioned) Christoffel transform of 
a cmc surface F : U C M 2 — > M 3 with unit normal field iV : U — > S 2 
is its parallel constant mean curvature surface F p = F + jjN. On the 
other hand, this parallel cmc surface is also a Darboux transform since a 
congruence of spheres with constant radius ^ is certainly a Ribaucour sphere 
congruence and — as is well known — a cmc surface F and its parallel cmc 
surface F p carry conformally equivalent metrics. In [|l^] we proved that 
this behaviour characterizes cmc surfaces: a surface F : U — > R 3 is cmc if 
and only if its (suitably scalecQ and positioned) Christoffel transform is a 
Darboux transform of F, too. 

Since we are working on Christoffel and Darboux transforms of isothermic 
surfaces it seems to be reasonable to use this characterization of cmc surfaces 
as a 

Definition (Discrete cmc net): An isothermic net F : T — > Ft 3 is called 
a discrete CMC net if a (suitably scaled^ and positioned) Christoffel trans- 
form F p : T — > R 3 of F is a Darboux transform of F , too. 

However, in Euclidean geometry this definition is rather unsatisfactory 
since it doesn't say anything about the mean curvature — which is claimed 
to be constant^: recall the meaning of "cmc". But, according to the above 
discussions we might define the (constant) mean curvature^ of a discrete 

13 Note that the sphere — which is certainly a surface of constant mean curvature — 
behaves exceptionally in this context: any Christoffel transform of the sphere is a minimal 
surface but its "parallel constant mean curvature surface" collapses to the center of the 
sphere. Therefore, for the sphere, "suitably scaled" means that the scaling factor j = 0. 

14 Because of the way we introduced the scaling factor A for the Christoffel transform 
in (0) this definition excludes the sphere: for the sphere the "suitably scaled" Christoffel 
transform would be the center of the sphere (cf. footnote |l3|), i.e. A = oo. 

In the geometry of similarities, however, this is a good definition: here, the value of the 
(constant) mean curvature is not an invariant. Note that the geometry of similarities arises 
much more naturally (by distinguishing a point at infinity) as a subgeometry of Mdbius 
geometry than the Euclidean geometry (where, additionally, a scaling of the conformal 
metric has to be choosen). 

16 In our discussions, Alexander Bobenko suggested a definition for the mean curvature 
function H : T — > R of an (arbitrary) discrete net F : T — > i? 3 at a point F mn of the 
net: consider the point C which has equal distances \F m _i^ n — C\ = \F m+ i^ n — C\ and 
|i*m,n-i — C| — \Fm,n+i—C\ to the pairs of opposite neighbours of F m>n and whose distance 
to F m , n satisfies \F m>n — C| 2 = ^(\F m+ i <n — C\ 2 + |F mj „ + i — C| 2 ). Its reciprocal distance 
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cmc net to be the reciprocal of the constant distance (vertex wise) of the 
nets F and F p — if we happen to prove the following 

Theorem: A Christoffel transform F c : T — > JR 3 of an isothermic net 
F : T — > R 3 is a Darboux transform of F , too, if and only if the distance 
\ F m, n ~ F m,n\ is constant. 

Since the quadrilaterals spanned by corresponding edgesQ of F and F c 
are trapezoids this theorem will be a consequence of 

The Trapezium lemma: The cross ratio of a trapezoid is real if and only 
if that trapezoid is isosceles. A quadrilateral (Qi, Q 2 , Q3, Q4) C H is an 
(isosceles) trapezoid, i.e. (Q1-Q4) \\ (Q2 — Q3) (and \Q\-Q 2 \ = \Q3-Q4D, 
if and only 2jQ 

DV(Qi, Q 2 , Q 3 , Qi) = -\ Qi -q^ Q \q^q 2 \2 ■ (23) 



It is clear that the four vertices of an isosceles trapezoid lie on a circle 
and consequently, its cross ratio has to be real. To understand the converse 
we assume that the vertices of a trapezoid lie on a circle. Then, the four 
vertices are obtained by intersecting this circle with two parallel lines. Thus, 
the trapezoid has a reflection symmetry showing that it has to be isosceles. 
This proves the first part of the trapezium lemma. Since a quadrilateral is 
uniquely determined by three of its points and its cross ratio the second part 
becomes clear by calculating the cross ratio (0) of an isosceles trapezoid. 

One direction in the proof of the theorem is now obvious: if the Christof- 
fel transform F c of an isothermic net F is also a Darboux transform, then, 

Hm,n = \F m ,n — C^ 1 can be considered the mean curvature of the net F at F m ^ n . Note, 
that C may not be a finite point — then the mean curvature vanishes (cf. the definition 
of discrete minimal surfaces in ||). 

17 Corresponding edges of F and any Christoffel transform F c are parallel — note that 
we restricted to the codimension 1 case. 

18 Note, that if \Qi — Q 2 \ < \Qi — Q3I then the trapezoid has legs Qi — Q2 and Q 3 — Q4 
and the isosceles trapezoid (Qi, Q2, Q3, Q4) is embedded. If \Qi — Q2I > \Qi — Q3 1 the 
edges Qi — Q2 and Q3 — Q4 become the diagonals of the trapezoid which is then not 
embedded. 
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all the trapezoids spanned by corresponding edges of F and F c are isosceles. 
Thus|3, \F mjn - F^ n \ = \F m+1}7l - = \F mjH+ i - F^ n+1 \. The other 

direction becomes clear, again, by calculating the cross ratio of an isosceles 
trapezoid: up to the "correct" scaling factor A c for the Christoffel trans- 
form F c , corresponding edges of F and F c have reciprocal lengthes @ and 
consequently, the cross ratios (0) become 



T)V(F F c F c F ) — \ C \F — F c I 2 

^ v \ A m,m 1 m>n s 1 m +\,ni 1 m+l,n) ^ \ 1 m,n 1 m ,n\ 

TN\Fm,ni F m n , F m n+ i, -F^n+l) A \F mn F m ^ n \ 



(24) 



— which proves the theorem. 

Now, we are able to give a definition of discrete cmc nets involving the 
mean curvature: 



Definition (Discrete net of constant mean curvature): An isothermic 
net F : V — > R 3 is a NET OF CONSTANT MEAN CURVATURE H if there is a 
Christoffel transform F p : T — > E? of F in constant distance^ 

\F nt m ~ F^ m \ = (25) 

F p is called the parallel cmc net of F. 

Note that this definition is symmetric in F and its Christoffel transform 
F p . So, by definition, the (correctly scaled) Christoffel transform of a discrete 
net of constant mean curvature if is a discrete net of constant mean curvature 
H, too. Clearly, rescaling of the Christoffel transform results in a discrete net 
of (generally) another constant mean curvature. Also note that the sign of 
the (constant) mean curvature H is not determined in the above definition^]. 

19 Moreover, after contemplating the relation between the lengthes of the edges and of 
the diagonals of an isosceles trapezoid we see that also the distances |F m „ — F^ +1 n \ and 
\F m . n — Ff n+l n \ are constant. 

20 As before (cf. footnote Q), the sphere is excluded by this definition. 

21 Choosing a sign for the constant mean curvature H would correspond to the choice 
of a unit normal field N : T — > S 2 . Note, that also the definition of the mean curvature 
function scetched in footnote |l6| includes no sign choice. 

From our previous discussions (cf. footnote [l9| ), it is clear that this mean curvature 
function (in the sense of footnote |l6|) , indeed, equals the constant H for a "discrete net of 
constant mean curvature H" (in the sense of the above definition). 
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From the work of Bianchi [Q] it is well known that (smooth) cmc surfaces 
allow oo 3 Darboux transforms in cmc surfaces. In ]12) we showed that a 
constant mean curvature Darboux transform F : U — > J? 3 of a constant 
mean curvature surface F : U C R 2 — > -R 3 has constant distance 

|F-Ff = i^^ (26) 

to its parallel constant mean curvature surface -F p . Herein, the mean cur- 
vature H c of F c of F is used to determine the scaling of the Christoffel 
transform in the Riccati equation (p|) . The 3-parameter family mentioned 
by Bianchi is obtained by solving the Riccati equation with an initial value 
satisfying (|26|). A similar theorem holds in the discrete case — we will obtain 
it as a consequence of the following 



Supplement to the Hexahedron lemma: //, in the hexahedron lemma, 
two of the faces (quadrilaterals) of the constructed hexahedron are (isosceles) 
trapezoids, then their opposite faces are (isosceles) trapezoids, too. 

To prove this supplement let — without loss of generality — the ini- 
tial quadrilateral {X\, X 2 , X 3 , X&) and the first constructed quadrilateral 
(Xi, X 2 , Z 2 ,Zi) be isosceles trapezoids^. Calculating the cross ratios of these 
two trapezoids yields 

n = + |Ma|2 



Xfj, 



Xi-x 4 \\X2-x 3 \ ^7) 

X\— Z1WX2— Zi\ 



|Xi-X 2 | 2 



- where all sign combinations can occur, depending on whether the two 
trapezoids are embedded or not. The remaining two points Z 3 and Z4 have 
to satisfy the condition 

DV(X U Z u Z 4 , X 4 ) = DV(X 2 , Z 2j Z 3 , X 3 ) = A. (28) 

In our situation, these two points can be constructed quite explicitely - 
showing that the quadrilaterals (X±, Zx, Z 4 , X4) and (X 2 , Z 2 , Z 3 , X 3 ) have 
parallel edges: without loss of generality we may assume that the circles de- 
termined by the point triples {X±, Xi,Zx} and {X 3 , X 2 , Z 2 } have the same 
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If two opposite quadrilaterals are trapezoids then there is nothing to prove. 
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Figure 8: A supplement to the hexahedron lemma 



radius. Since (X4 — Xi) \\ (X 3 — X 2 ) and (Xi — Z%) \\ (X 2 — Z 2 ) we con- 
clude IX4 — Zi\ = \X 3 — Z 2 \ - this argument is scetched in figure [|, in 
case the first two trapezoids are embedded and in case the initial one is 
embedded and the second one is not. Thus, we find a reflection which in- 
terchanges X4 with Z 2 and X3 with Z\. Now, the two points we had to 
construct are the images Z4 of X 2 and Z 3 of X\ — it is easily checked that 
the quadrilaterals (Xi, Z±, Z4, X 4 ) and (X 2 , Z 2 , Z 3 , X 3 ) have the correct cross 
ratio. Since these quadrilaterals obviously have parallel edges all quadrilat- 
erals (X u X 2 ,X 3 ,Xi), (Xi,X 2 ,Z 2 ,Zi), (Zx, Z 2 , Z 3 , Z 4 ) and (X 4 , X 3 , Z 3 , Z 4 ) 
are trapezoids and since their cross ratios are real they are isosceles. This 
proves the above supplement to the hexahedron lemma. 
With this knowledge we can now prove the announced 



Theorem: If F : T —> ImH is a discrete net of constant mean curvature 
H , then any solution F : Y — > ImH of the Riccati type system ftTBj ) with the 
initial condition 

\F mo ,n — Fm ,n \ = Jp ~ J^j (29) 

is a discrete net of constant mean curvature H . Herein, X p is the parameter 
of the Darboux transform which transforms the cmc net F into its parallel 
cmc net F p . 

To prove this theorem we will have to construct the parallel cmc net F p of 
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F. This net is a A p -Darboux transform of F and a A-Darboux transform of the 
parallel cmc net F p of F at the same time — consequently, the construction 
will be similar to that in the proof of the permutability theorem for the 
Darboux transform (cf. Fig. ^|). 

By the trapezium lemma, the quadrilateral (F moM ,^ OM , F m ,n i F m ,n ) 
which is first constructed is an isosceles trapezoid because of (f29"|). Since all 
the quadrilaterals spanned by corresponding edges of F and its parallel cmc 
net F p are also isosceles trapezoids, according to our supplement to the 
hexahedron lemma, there is a net F p : T — > ImH whose edges are parallel to 
the corresponding edges of F and which has constant distance r^r to F. The 

only thing left is to show that F p is a Christoffel transform of F — which 
becomes clear by contemplating the cross ratios 



DV(-F mjTl , F^ n , F m+ln , -F m +i iTC ) A p 



(30) 



since |F m>n - F^J 2 = ^2 it follows 



{Fm+l,n F m,n){,F m +l,n -^m,n) \PH 2 

(F m<n+ i — F m<n )( K F m n+1 — -F^ n ) — — \PH 2 ' 



(31) 



Thus, F p is a Christoffel transform of F and consequently it is its parallel 
cmc net in distance rjU. This completes the proof. 

Now, let's assume we have two cmc Darboux transforms F\ and F2 of a 
discrete net F of constant mean curvature H, and their parallel cmc nets 
F p , Ff and F% which are A p -Darboux transforms of the original nets. By the 
"fancy" version of Bianchi's permutability theorem, the picture gets com- 
pleted with two discrete isothermic nets F and F p — F being a A2-Darboux 
transform of Fx and a Ai-Darboux transform of F 2 and F p being a A p -Darboux 
transform of F, a Ai-Darboux transform of and a A2-Darboux transform 
of Ff. Since, moreover, F p and Ff are the parallel cmc nets of F and 
Fi the quadrilaterals (F, F p , Ff, Fi) and (F, F p , F%, F 2 ) are isosceles trape- 
zoids. Thus, by our supplement to the hexahedron lemma, the quadrilater- 
als (Fi, Ff, F p , F) and (F 2 , F| ', F p , F) are isosceles trapezoids, too. Conse- 
quently, the net F p is the parallel cmc net of the cmc net F — we just proved 
the following permutability theorem for cmc Darboux transforms of discrete 
cmc nets (cf. M): 
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Theorem: If F\,F 2 : T — > ImH are two cmc Darboux transforms of a 
discrete net F : Y — > ImH of constant mean curvature H with parameters 
Ai and A 2 , respectively, then there exists a net F : T — > Iml/ 0/ constant 
mean curvature H which is a X 2 -Darboux transform of Fi and a Xi-Darboux 
transform of F 2 at the same time. The nets F, F 1 , F 2 and F have constant 
cross ratio ^ = DV(F, F 2 , F, A)- 
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